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1. Introduction and physics 

Neutrinos play a crucial role in modern particle, nuclear and astrophysics including cos- 
mology. It has been the major achievement of the last 15 years to show that neutrinos 
have a non-vanishing rest mass. The evidence arises from a deficit of upward going at- 
mospheric muon neutrinos confirmed by long baseline accelerator experiments and the 
solution of the solar neutrino problem being confirmed by nuclear reactors. All can be ex- 
plained by neutrino oscillations, which are depending mass differences Amfj = m 2 — mf 
assuming a two flavour mixing only. The determination of absolute neutrino masses is 
now a major issue, because neutrino oscillation experiments do not allow this. The classi- 
cal way search for a rest mass of the neutrino is the study of the endpoint region of electron 
spectra in beta decay (see [ 1 1 for a recent review). The KATRIN experiment is well on 
its way to improve the current bound of about 2.2 eV for D e by an order of magnitude. 
Further bounds on the total sum of neutrino masses can be obtained from cosmological 
studies. Another laboratory process is the rare nuclear decay of neutrinoless double beta 
decay. 

(Z,A) -> (Z + 2,A) + 2e- (Oz^/3/3-decay ) (1) 

Single beta decay must be forbidden or at least strongly suppressed to observe this decay 
and 35 potential isotopes exist in nature. As can be seen from Equation[T]the given decay 
mode is violating total lepton number by two units and thus not allowed in the Standard 
Model. Being a decay the observable is a half-life which can be linked to the quantity of 
interest e via 

(T^)- 1 = G PS I M Nuc | 2 e 2 (2) 

with Gps being the phase space, | Mn uc | the involved nuclear matrix element for the 
physics process considered to describe this decay and e the quantity of interest. Various 
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BSM processes can be considered among them light and heavy Majorana neutrino ex- 
change , right-handed weak currents, R-parity violating SUSY (A' m ) and double charged 
Higgs bosons. If O^/3/3-decay is ever observed it was shown that this would imply that 
neutrinos are Majorana particles 0. However, given all the possible processes it will 
become an important question what the individual contributions of the considered pro- 
cesses will be. The LHC will help to restrict this by performing searches for new parti- 
cles in the TeV range. An example how this complementary information from LHC and 
O^/3/3-decay can be used within the context of left-right symmetric theories is given in 
. For a recent extensive review on the particle physics in double decay see H) . 
The standard interpretation considered here is the one using light Majorana neutrino ex- 
change. In this case e is the effective Majorana neutrino mass (m ee ) given by 

£ = (m ee ) =| ^2 U ei m i I (3) 
i 

with as the mixing matrix elements containing the electron neutrino. 
2. Double beta decay and neutrino oscillation results 

For the following discussion a restriction to the light Majorana neutrino case is done. It 
is evident from Eq.|3]that the expectation for (m ee ) in double beta decay depends on the 
neutrino oscillation parameters. It should be noted that the mixing matrix of relevance is 
given by 

U = U PMN s(0 12 ,0 13 ,9 23 ,e lS ) x diag(l,e ia ,e i0 ) (4) 

with the standard leptonic mixing matrix Upmns ar, d two additional CP-phases a and 
(3, called Majorana phases, which appear if neutrinos are their own antiparticle. These 
phases do not show up in oscillation experiments. Hence, (m ee ) can be written as a sum 
of three terms 

(me.) =| ml I + I m 2 e | e 2ia + \ ml \ e 2lfi (5) 
with the individual contributions given as 

m\ =| U e \ | 2 mi — mi cos 2 9\ 2 cos 2 #13 

m 2 e =| U e2 | 2 m 2 = m 2 sin 2 9 12 cos 2 6>i 3 (6) 
ml =| U e3 | 2 m 3 = m 3 sin 2 6 13 

Latest global fits to available oscillation parameters are given in |5]|6|. An important 
new ingredient is that 3-flavour fits to solar neutrinos [7|, observations from T2K (8) and 
MINOS (9) as well as first Double Chooz results [ 10 1 favour a non-zero value of 6*13. As 
the oscillations do not fit the absolute scale there are two options for arranging the mass 
eigenstates, either m 3 > m 2 > mi (normal hierarchy, NH) or m 2 > mi > m 3 (inverse 
hierarchyJH). If the neutrino masses turn out to be close the the current limit from beta 
decay there is a quasidegeneracy (mi « m 2 ~ m 3 = m ). In case of hierarchies the two 
larger masses can be expressed as 



m 2 = %J m\ + Am,Q m 3 = m\ + Am 2 trra (normal) 
m 2 = Jm% + Am| + Am 2 tm m x = yjm% + Am 2 tm (inverted) 
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(7) 

with the solar splitting Am| = 7.59^Q 2 g x 10~ b eV 2 and the atmospheric splitting 
Am 2 atm = 2.49 ± 0.09 x lO^eF 2 (normal) , Am 2 atm = — 2.343±g;Jg x 10- 3 eF 2 (in- 
verted) 1 5 1. The principal behaviour of (m ee ) as a function of the lightest mass eigenstate 
is shown in Figure [l]for an exemplaric value of sin 2 26*13 = 0.02 and the dependence of 
individual parts on the parameters. For (m ee ) values larger than about 100 meV neu- 




Figure 1. The effective Majorana neutrino mass (m ee ) as a function of the smallest 
neutrino mass and the principal dependence of individual features on the oscillation 
parameters. A value of sin 2 28n = 0.02 has been assumed, furthermore i? 2 = tan 2 #12. 
Uncertainties introduced due to the involved nuclear matrix elements are not included 
(from (TTJ). 



trinos are almost degenerate, the inverted hierarchy covers a range between about 10-50 
meV and below 10 meV is the region of the normal hierachy. As can be seen, in the NH 
there can be a cancellation among the terms, the allowed region for that becomes larger 
with increasing sin 2 26*13. There is no such effect in the IH because of the non-maximal 
solar mixing angle 6*12- In addition the gap between the two bands in the hierarchical re- 
gion will shrink with increasing sin 2 26*13. Half-lives for the IH are in the region beyond 
10 26 years while half-lives in the NH are well beyond 10 28 years. 



3. General experimental considerations 

Evidently measurements of half-lives well beyond 10 25 years are by now means trivial. 
As signal for the process given in Equation ([T]i serves a peak in the sum energy spectrum 
of the two electrons equivalent to the Q-value of the nuclear transition. The corresponding 
half-life in case of no background is given by the radioactive decay law 

Tl" = In 2matN A /Npp (8) 

with m the used mass, a the isotopical abundance of the double beta emitter, t the mea- 
suring time, Na the Avogadro constant and Npp the number of double beta events, which 
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has to be taken from the experiment. If no peak is observed and a constant background 
(i.e. all potential energy depositions in the region of interest, i.e. around the Q-value, 
not being neutrinoless double beta decay) is assumed scaling linearly with time, Npp is 
derived as 



1 l/2 



M x t 
B x AE 



(9) 



where e is the efficiency for detection of the total energy of both electrons, AE is the 
energy resolution at the peak position and B the background index normally given in 
counts/keV/kg/year. Hence, the most crucial parameters are a high detection efficiency 
and high abundance of the isotope of interest. This is the reason why almost all next 
generation experiments are using enriched materials and the "source = detector" approach. 
Furthermore, the energy resolutiorQ should be as good as possible to concentrate the few 
expected events in a small region and ideally the experiment should be background free. 
An irreducible background is the Standard Model process 2^/3/3-decay 

(Z,A) -> (Z + 2,A) + 2er +2v e (2;//3/3-decay ). (10) 

Here again energy resolution matters, because of the continuous spectrum of the 2v (3j3 '-decay mode 
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Figure 2. Schematic plot of the sum energy spectrum of the two electrons in double 
beta decay. O^/3/3-decay results in a peak at the Q-value of the transition. Various 
modes can be characterised by the phase space dependence (Q-E) n . The mode n=5 is 
the 2^/3/3-decay while the modes n=l,3,7 involve the emission of a majoron, a Gold- 
stone boson linked to the spontaneous breaking of lepton number. The different modes 
belong to different behaviours of the majoron with respect to weak isospin. The indi- 
vidual contributions are not to scale. 



it is only its high energy part leaking into the peak region (see Figure |2])). Nevertheless, 

'Care must be taken as for traditional reasons different detector technologies use either the Gaussian a 
or the Full Width at Half Maximum AE to quote energy resolution. The relation among the quantities is 
AE = 2.35a. 



4 



Kai Zuber 



this can be a worry as the half-life is typically several orders of magnitude shorter than the 
expected one for O^/3/3-decay . As the decay rate for Oz^/3/3-decay scales with Q 5 only 
isotopes with Q-values above 2 MeV are considered for experimental searches. They are 
listed together with some important numbers in Table ([TJ. 



Isotope 


nat. abund. 

(%) 


Q-value 
(keV) 


rp'lV 

1 l/2 

(10 20 yrs) 


Experiment 


48 Ca 


0.187 


4272 ± 4 


44+U-ue 


CANDLES 


76 Ge 


7.8 


2039.006 ± 0.050 


15 ± 1 


GERDA, MAJORANA 


82 Se 


9.2 


2995.5 ±1.9 


0.92 ± 0.07 


SuperNEMO, LUCIFER 


96 Zr 


2.8 


3347.7 ± 2.2 


0.23 ± 0.02 




100 Mo 


9.6 


3034.40 ±0.17 


0.071 ±0.004 


AMoRE 


no pd 


11.8 


2017.85 ± 0.64 






116 Cd 


7.5 


2813.50 ±0.13 


0.28 ± 0.02 


COBRA, CdW0 4 


124 Sn 


5.64 


2287.8 ± 1.5 






130 Te 


34.5 


2527.518 ±0.013 


6.8 t\i 


CUORE 


136 Xe 


8.9 


2457.83 ± 0.37 


21.1 ±2.5 


EXO, NEXT 










KamLAND-Zen 


i50 Nd 


5.6 


3371.38 ±0.20 


0.082 ± 0.009 


SNO+, DCBA 



Table 1. Table showing the eleven candidate isotopes with a Q-value larger than 2 
MeV. Given are the natural abundances, Q-values as determined from precise Penning 
trap measurements (those with sub keV errors) or from the Atomic Mass Evaluation 
2003 [ 12 1, the measured averaged 2i//3f3 -decay half lives as recommended in [ 13] plus 
the recent measurement of 136 Xe 1141 . The last column shows the experiments ad- 
dressing the measurement of the corresponding isotope. For some experiments only 
the "default" isotope is mentioned as they have the option of exploring several ones. 
Several additional research and development projects are ongoing. 

As mentioned, most experiments follow the approach that the source is equal to the 
detector, i.e. building a detector which contains the isotope of interest. Technologies 
used for that are semiconductors, cryogenic bolometers, scintillators and liquid noble gas 
detectors. The alternative is to use tracking devices in form of TPCs containing thin foils 
of double beta emitters. Here single electron spectra and opening angles can be measured 
as well. 

4. Experimental status 

Three different goals are considered for future investigation depending on the outcome of 
the individual steps. The first one is to probe the claimed observation of O^/3/3-decay in 
76 Ge with a half-life of 2. 23±0.04xl0 25 yrs EH El. If this is confirmed, the next genera- 
tion will collect sufficient statistics for a precision half-life measurement. It might be even 
possible to perform an intrinsic consistency check by also looking at the O^/3/3-decay de- 
cay into the first excited state 1 18 1. Furthermore, due to the "multi-isotope" approach an 
observation in one isotope predicts a half-life for the other ones by taking into account 
the uncertainties in nuclear matrix elements, see section ([6]) . Observing peaks at different 
positions in energies within the right range of half-lives excludes potential unrecognised 
backgrounds. Anyhow, the "multi-isotope" ansatz is needed to compensate for matrix 
element uncertainties. 

In case the evidence is not confirmed, the next goal must be the region of the inverted 
hierarchy, i.e (m ee ) below w 50 meV. The neccessary half-life requirement to touch this 
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region is given in |[T9l . As the requirements on mass and background for this purpose 
are already demanding, fully excluding the IH is more than challenging. If there is no 
signal found in the IH the final is the exploration of the normal hierarchy. However, 
for this ton scale experiments with extraordinary low background have to be considered 
and completely new background components have to taken into account, for example 
neutrino-electron scattering due to solar neutrinos lf20l . For that half-lives well beyond 
10 28 yrs have to be measured. 

Currently the field is in the transition towards the next generation of experiments. The 
year 201 1 is seeing the start of four new double beta experiments, namely GERDA (using 
76 Ge ), EXO and KamLAND-Zen (using 136 Xe ) and CANDLES (using 48 Ca ). GERDA 
ETIl is a next generation experiment based on Ge-semiconductors containing 76 Ge based 
in the Gran Sasso Underground Laboratory (Italy). The idea is to run bare crystals within 
LAr which serves as shielding and cooling. In a first phase a total of eight isotopically 
enriched detectors (from the former Heidelberg-Moscow and IGEX experiments) with 
a total mass of 17.7 kg have been deployed and official data taking started on 1. Nov. 
201 1. A spectrum from a commissioning run using 3 enriched detectors is shown in Fig- 
ure (|3j. As can be seen the 2i//?/3-decay contribution is already clearly visible. Currently, 
for a second phase another 35.5 kg of enriched Ge is in the process of conversion into 
BEGe detectors. These point-contact detectors allow a optimised separation of single 
site energy depositions (like double beta decay) and multiple site interactions (most of 
the background like Compton scattering with an additional second interaction) by using 
pulse shape analysis |22"I . 

Due to the relative cheap enrichment of noble gases, two large scale experiments based 
on enriched 136 Xe have started. First is EXO-200, using about 175 kg LXe in form of a 
TPC located at WIPP (USA). The early data look very promising (Figure [3} and a first 
half-life for the 2^/3/3-decay of 136 Xe could be determined, see Table (JT|). A unique op- 
tion explored for the future would be the detection of the daughter ion which should result 
in major background reduction. A second Xe approach is KamLAND-Zen using the well 
understood infrastructure of the KamLAND experiment. To perform double beta decay 
with Xe-loaded liquid scintillator a special mini-balloon was constructed and deployed 
within KamLAND. In this way about 330 kg of enriched Xe could be filled in the detec- 
tor and data taking has started in September 2011. 

Last but not least there is CANDLES using 305 kg of CaF2 scintillators, focussing on 
48 Ca , the isotope with the highest Q-value of all double beta emitters. The experiment is 
installed in the Kamioka mine (Japan) and data taking has started recently. 

5. Future experiments and ideas 

With the experiments mentioned in the previous session running, there will be more online 
in a few years time and additionally still interesting research and development is ongoing, 
a few of them are mentioned here, see Table ([TJ. From the semiconductor point of view, 
two more projects exist, MAJORANA and COBRA. MAJORANA is using Ge-detectors 
as GERDA and is considering to have up to 30 kg enriched detectors running by 2014 
in the Sanford Underground Laboratory (USA). COBRA is exploring CdZnTe semicon- 
ductor for the search of 116 Cd . The detectors are much smaller than the typically used 
Ge-diodes thus having a much higher granularity and discrimination against multi-site 
events will be realised by multi-detector events. Furthermore, for practical purposes it 
is convenient that these detectors will run at room temperature. A unique feature which 
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Low energy spectrum with enriched HPGe detectors 




Energy (keV) reconstructed energy {5|i <keV> 

Figure 3. Left: Spectrum from a commissioning run of GERDA using 3 enriched 
Ge-detectors. The three major visible components are the 39 Ar - decay, a 42 K contribu- 
tion due to 42 Ar decay and the 2v/3f3-decay of 76 Ge . Right: Initial results of EXO-200 
lead to a first discovery and half-life measurement of the 2^/3/3-decay of 136 Xe being 
the dominant contribution to the spectrum (from 1 15], details can be found in 1 14 1). 



is explored is pixelisation which should allow massive background reduction by particle 
identification due to tracking. A prototype detector of about 0.5 kg is currently installed 
in the Gran Sasso Laboratory. 

In a similar spirit as KamLAND-Zen another large scale scintillator experiment is in the 
building up phase. SNO+ aims to use in a second phase of operation (a first phase should 
be pure scintillator mainly for the study of solar neutrinos) 1000 1 of Nd-loaded scintillator 
for the search of 150 Nd decay. The default assumption is a 0.1% loading of the scintil- 
lator, resulting in an amount of 43 kg of 150 Nd , while an optimisation of 0.3% loading 
is under investigation. 150 Nd has the second highest Q-value of all double beta isotopes 
(3371 keV) and was always among the theoretically most preferred isotopes. SNO+ is 
supposed to start in 2013. Other groups explore various solid state scintillators, for exam- 
ple CdW0 4 . 

A different class of detectors not mentioned yet but with a lot of options are cryogenic 
bolometers. The most advanced approach is CUORE using 750 kg of TeC>2 crystals to 
search for the 130 Te decay. The benefit of using Te is the high natural abundance of 
130 Te . This experiment is in the building up phase at Gran Sasso Laboratory and is 
planned to start data taking with the full amount in 2014. Various other bolometers are 
currently explored like CaMo0 4 (AMoRE) and ZnSe (LUCIFER) to search for 100 Mo 
and 82 Se respectively. 

Finally tracking devices in form of TPCs using the double beta emitter in form of thin 
foils are explored. Due to the tracking it is possible to measure the individual electron 
energies and the opening angle between them. This might be important if Of /3/3-decay is 
discovered as the neutrino mass mechanism and right-handed weak current contribution 
differ significantly in these quantities. As a next step in a serious of experiments Su- 
perNEMO is planning to use at least 100 kg of 82 Se in form of 20 TPC modules with 5 
kg source mass per module. A first demonstrator module is supposed to start data taking 
in 2014. 

6. Nuclear matrix elements - Theory and experiment 

The conversion of an observed half-life or its limit into (m ee ) requires the knowledge of 
the nuclear transition matrix elements, see Equation Any uncertainty in the nuclear 
matrix element due to its square dependence results in a significant uncertainty in the de- 
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duced neutrino mass. This is another important argument for the multi-isotope approach. 
Various theoretical methods have been applied in the past, dominantly nuclear shell model 
(NSM) and quasi-random phase approximation (QRPA) calculations, which were recently 
joined by new methods like the Interaction Boson Model (IBM) and energy density func- 
tional treatment (GCM). A compilation of the available calculations normalised to one 
value of ro and gA is given in Figure Q. As can be seen there is quite some spread in the 
obtained values and uncertainties are not claimed for all methods. 

As all double beta emitters are even-even nuclei their ground state transitions are char- 
acterised as + — > + transitions. Given that fact that almost all 2^/3/3-decay half-lives 
are known the associated matrix element M 2v can be deduced. However, these are pure 
Gamow-Teller transitions only mediated by the 1 + states of the intermediate nucleus and 
restricted the less than 5 MeV and they contribute only a fraction to M 0lJ . In the latter 
case all levels up to about 100 MeV with all multipolarities can contribute. Furthermore, 
effects and treatment of short range correlations become very important l23l . The matrix 
elements itself must be calculated separately for the individual processes discussed in sec- 
tion ([TJ. Thus, they contain information on the physics process involved as well ll24l|26l . 
Given the complexity of the problem and a lack of missing experimental information, a 
program was started to improve the situation from the experimental point ll27l . Part of 
the program is precise Q-value determinations using Penning traps (see Table [TJ, charge 
exchange reactions in form of ( 3 He, t) and (d, 2 He) reactions to determine the Gamow- 
Teller strength Bqt for the transitions to the intermediate l + -states and nucleon transfer 
reactions. This lead to new insights and refinements of the calculations. 
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Figure 4. Compilation of nuclear matrix element calculations for the eleven relevant 
isotopes. Existing publications were used and normalised to a nucleon radius of ro = 
1.2 fm and £a=1.25 (from fl9l ). 



8 



Kai Zuber 



7. Alternative processes including positrons and electron capture 



An equivalent process to the one discussed is /3 + /3 + -decay also in combination with elec- 
tron capture (EC). There are three different variants possible depending on the Q-value: 



2e~, 



(Z,A) 
(Z,A) 
(Z,A) 



(Z-2,A)+2e+(+2v e ) (/3+f3+) 
(Z-2,A) + e+(+2v e ) (13+ /EC) 
(Z -2,A)(+2u e ) (EC /EC) 



(11) 
(12) 

(13) 



f3 + f3 + is always accompanied by EC/EC or /3 + /EC-decay. The positron production 
reduces the effective Q-value by 2m e c 2 per positron. Therefore, the rate for /3 + /3 + is 
small and energetically only possible for six nuclides, however it would have a striking 
signature with four 511 keV gamma rays. It was shown that the /3 + /EC-mode has an 
enhanced sensitivity to right-handed weak currents |28 1 and might be valuable to explore 
if O^/3/3-decay is discovered. The full Q-value is available in the EC /EC mode which is 
the hardest to detect experimentally. However, it was proposed [29 30 1 that if an excited 
state of the daughter nucleus is degenerate with the original ground state a resonance 
enhancement in the decay rate could occur and the de-excitation gammas would serve as 
a nice signal. Due to the sharpness of the resonance a more detailed study of candidates 
had to wait for Penning traps entering the field and exploring reasonable candidates. The 
most reliable one seems to be 152 Gd (see Figure|5]l where such a scenario is realised BP . 
Despite this nice effect, to achieve the same sensitivity of (m ee ) as in O^/3/3-decay seems 
to require a measurement an order of magnitude longer making this method slightly less 
attractive. 




™ " 10 2 
o £ 
d) ro 



Q -B, h / keV 

Figure 5. Resonance enhancement factor for neutrinoless ECEC relative to 54 Fe based 
on precision mass spectrometry with Penning traps (from [31 1). 
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